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Abstract 

We consider the conformal Einstein equations for 1 < 7 < 2 poly- 
tropic perfect fluid cosmologies which admit an isotropic singularity. 
For 1 < 7 < 2 it is shown that the Cauchy problem for these equations 
is well-posed, that is to say that solutions exist, are unique and depend 
smoothly on the data, with data consisting of simply the 3-metric of 
the singularity. The analogous result for 7=1 (dust) is obtained when 
Bianchi type symmetry is assumed. 



1 Introduction 

Isotropic singularities are a class of cosmological singularities which are both 
physically interesting and mathematically tractable (Goode and Wainwright 
1985, Tod 1987, Tod 1992). Put simply, an isotropic singularity is one which 
can be removed by rescaling the spacetime metric with a single function 
which becomes singular on a spacelike hypersurface. 

The motivation for studying these singularities comes from Penrose's 
Weyl Tensor Hypothesis (WTH) (Penrose 1979, Penrose 1990). Penrose 
hypothesised that initial singularities, notably the Big Bang, should be dif- 
ferent in character from final singularities, as formed in gravitational collapse 
or the Big Crunch. Specifically (Penrose 1979), he suggested: 

I propose that there should be a complete lack of chaos in the 
initial geometry. . . This restriction on the early geometry should 
be something like: the Weyl curvature Cabcd vanishes at any 
initial singularity. 
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Penrose expected this condition to imply that the subsequent evolution 
would necessarily be close to a Friedman-Robertson- Walker (FRW) model. 
He was led to the above hypothesis by the need for some kind of low entropy 
constraint on the initial state of the universe, at the same time as the matter 
content was in thermal equilibrium. He argued firstly that the need for a 
low entropy constraint follows from the existence of a second law of thermo- 
dynamics and secondly that low entropy in the gravitational field must be 
tied to constraints on the Weyl curvature. Penrose went on to suggest that 
this constraint on the initial Weyl tensor should be a consequence of an as 
yet undiscovered quantum theory of gravity. 

In order to make mathematical progress with the WTH one must first 
say exactly what is meant by a 'finite Weyl tensor' singularity, a singularity 
at which the Weyl tensor is finite while, of necessity, the Ricci tensor is 
not. There may be more than one plausible way to do this. One strategy, 
and the one we shall adopt, is first to make the following definition: A 
spacetime {M,gab) is said to admit an isotropic singularity if there exists a 
manifold with boundary M D M, & regular Lorentz metric gab on M, and a 
function CI defined on M, such that 

9ab = ^^9ab for n>0 (1) 

17 ^ on S (2) 

where S is a smooth spacclikc hypcrsurfacc in M, called the singularity sur- 
face. Since the Weyl tensor with its indices arranged as C"{,cd conformally 
invariant, isotropic singularities form a well-defined class of cosmological 
singularities with finite Weyl tensor. 

Consideration of the conformal Einstein equations near an isotropic sin- 
gularity, for various matter models, leads naturally to a class of singular 
Cauchy problems for the unphysical metric gab, with data given on S, for 
which one may seek to prove suitable existence and uniqueness theorems. 
This program was begun in (Tod 1990, 1991), where the problem was treated 
in the case of a perfect fiuid with polytropic equation of state as matter 
source. There the conformal field equations were written as an evolution 
system for gab with data just the 3-metric of the singularity surface. It was 
noted in particular that a uniqueness result for these equations would prove 
an earlier 'Vanishing Weyl Tensor' conjecture, (Tod 1987), which said that, 
given an equation of state, the Weyl tensor must vanish everywhere in M if 
it vanishes at an initial isotropic singularity. 

Newman (1993a, b) considered cosmologies with an isentropic perfect 
fiuid as source and for which the conformal factor $7 remained a smooth 
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coordinate at S. This smoothness condition imphed that the equation of 
state approached that for radiation (7 = 4/3) at the singularity. By tak- 
ing 7 = 4/3 throughout M Newman was able, with a careful choice of field 
variables, to write the equations in symmetric hyperbolic form but with a 
singular 'j' forcing term. Again the free data at the singularity consisted 
of just the 3-metric there. As there was no existence and uniqueness theo- 
rem available for the singular evolution equations obtained, a special study 
was undertaken by Newman, the result appearing in (Claudel and Newman 
1998). We shall refer to the existence and uniqueness theorem obtained 
there as the Newman-Claudel Theorem. The conclusion is that there is ex- 
actly one 7 = 4/3 perfect fluid cosmology with an isotropic singularity for 
each 3-metric given on S. 

Our aim in this paper is to show that by imposing a suitable 7-dependent 
differentiability condition on the conformal factor 0,, a similar analysis to 
that of Newman can be performed for polytropic perfect fluids with index 
in the range 1 < 7 < 2. It will then follow from the Newman-Claudel Theo- 
rem that we have existence and uniqueness of a perfect fluid cosmology for 
each 7 in this range, for each 3-metric on E. 

The plan of the paper is as follows: 

In section 2 we set down our conventions and give a few notes on con- 
formal rescaling. 

An account of the theory of singTilar symmetric hyperbolic systems de- 
veloped by Claudel and Newman (1998) is given in section 3. 

In section 4 we consider cosmologies with 1 < 7 < 2. Using comoving 
coordinates as in (Newman 1993a,b) we are able to define a set of variables 
in terms of which the conformal field equations can be written in a suitable 
singular symmetric hyperbolic form. The results of (Claudel and Newman 
1998) can be applied and the Newman-Claudel theorem then gives existence 
and uniqueness of solution. 

The case 7 = 2 has to be treated separately and this we do in section 
5. We use harmonic rather than comoving coordinates and this enables us 
once again to write the field equations in the form required in (Claudel and 
Newman 1998) and to conclude existence and uniqueness of solutions from 
the Newman-Claudel theorem. 

Finally in section 6 we consider the case 7 = 1. Here it has not been 
possible to write the conformal equations in full generality in symmetric 
hyperbolic form and for this reason we are forced to restrict attention to 
spatially-homogeneous solutions with Bianchi type symmetry. Subject to 
this restriction, we are able to prove existence and uniqueness of solutions. 
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given an initial, spatially-homogeneous 3-metric, by an extension of a theo- 
rem from (Rendall and Schmidt 1991). 



2 Conventions and conformal rescaling 

Throughout the paper we take the spacetime metric to have signature 

(+-—)■ 

For a metric gab our definition of the Riemann tensor is 

iyaVb-VjNa)V' = R''^by^ (3) 

while for the Ricci tensor we take 

Rab = R^acb (4) 

Now whenever we're considering a rescaling as in (l)-(2), tilded quantities 
will always refer to the singular, physical spacetime (M, gab), while un-tilded 
quantities will refer to the regular, unphysical spacetime {M,gab)- For met- 
rics gab and gab related by (1), the Ricci tensors are related by: 

Rab = Rab - 2V„ Vb log n + 2V„ log nVb log n 

- 5a6 (n log n + 2Vc log logQ) (5) 

and the Weyl tensors by: 

Cfocd = C^cd (6) 

If i"' is a unit vector in M then = is one in M, and if Wa is a unit 
covector in M, then Wa = ^~^Wa is one in M. With respect to a choice of 
unit timelike vector t"", the electric and magnetic parts of the Weyl tensor 
of gab are defined, respectively, by: 

Kb = C'adbict'^ Bab = 2^ae'^C{^tHf (7) 

The projection hab orthogonal to is defined by 

hab = 9ab - ^aU (8) 

If is hypersurface orthogonal one then calculates that 

Eab = -^Rab + KKab " KaeK% + ^h^h'Rsr 
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+ (^Ih^'^Rcd - Kb (9) 

where ^Rab is the Ricci tensor of hab, Kab is the second fundamental form 

of the surfaces orthogonal to and K = h'^'^Kab- 

Also 

B,b = DaK,if,e^f (10) 

where Da is the intrinsic covariant derivative associated with hab a^nd 6^60 is 
the volume element associated with hab- 
By (6)- (7) we have 

Eab = Eab Bab = Bab (H) 

It follows from (1) that, if Vt vanishes at a surface S in M then S is a 
singularity in M but, from (11), Eab and Bab are finite at the singularity. 
It follows from (10) that if S has vanishing second fundamental form Kab 
then Bab is automatically zero at S and from (9), subject to the term in the 
4-dimensional Ricci tensor being a multiple of hab, that Eab also vanishes at 
S if and only if S is a 3-dimensional space of constant curvature. This will 
be significant in section 5.5. 

3 Singular evolution equations 

The fundamental technique used by Newman (1993b) in tackling the confor- 
mal Cauchy problem for the Einstein-perfect fluid equations was to choose 
certain field variables so that the resulting equations appeared in symmetric 
hyperbolic form (as in equation (15) below). Although the system was hy- 
perbolic, the unavoidable ' j' singularity in these equations precluded the use 
of standard existence and uniqueness theorems as found in e.g. (Kato 1975, 
Rackc 1992). A special study was undertaken by Newman to deal with the 
singularity, and he was able to prove well-posedness of the CaTichy problem 
for a class of singular evolution equations into which the 7 = 4/3 equations 
fell (Claudel and Newman 1998). 

We will show in section 4 that the conformal field equations for other 7 in 
the range 1 < 7 < 2 can, after a little algebra, be forced into the class studied 
in (Claudel and Newman 1998), giving us the result we are looking for. 

For ease of reference we now present the main theorem of (Claudel and 
Newman 1998), both in its abstract form and its application to symmetric 
hyperbolic systems. 



5 



Theorem 1 (Newman-Claudel 1998). Let 

^ . . . ^ X° (12) 

be a sequence of continuous dense inclusions of reflexive Banach spaces, and 
let X = X°. 

Consider the quasi-linear, first order evolution problem 

^tM^ut = At{ut)ut + Mut) + ^Ft{ut)ut, 0<t<T, 

u{0) = u° (13) 

in an open neighbourhood U of € X for families {ftix) G X}{t,x)£[o,T]xU-: 
{Ft{x) G B(X)}(j 2:)e[o,r]x[/- i^i^) being the space of bounded operators 
in X), and a family {^t(a:)}(t^2:)G[o,T]xC/ of operators in X such that 

(i) XP+i C D{{At{x))P) for all (i,x) G [0,r] x U; 

(ii) n^Gn-l^'^^; 

(iii) Ft(x)^x° = for all {t,x) G [0,T] x ?7. 

Let C U he an open neighbourhood of in and, for each p > 
0, let C/P+^ be an open neighbourhood of contained in U^. Suppose 
there exists /? > and, for each (t, x) G [0, T] x U, an equivalent renorm- 
ing X(j.^^^ of X, such that 

(iv) At(y) G G(X(j,j,),l,/3) for aU (t,y) G [0,r] x (i.e At{y) gener- 
ates a contraction Co semigroup on ) ; 

(v) [0,T] X U 3 X 1-^ \\ llx(t^) is Lipschitz continuous; 

(vi) for each t G [0, T] and allp > 0, UP 3 x ^ ft{x) is in Lip{UP, XP), UP 3 
X ^ Ft{x)\XP is in Lip{UP,B{XP)) and UP 3 x ^ At{x)\XP+^ is in 
Lip{UP, B{XP'^^,XP)), with f-independent Lipschitz constants in each case; 

(vii) for all p, r : < r < p one has /|[0,r] x UP'"' G C''{[0,T] X 
UP-'',XP^''), and that [0,r] x UP-'' 3 {t,x) ^ Ft{x)\XP-'' is in C''{[0,T] x 
C/^'-^5(X^'-0), and that [0,r] x UP-'' 3 {t,x) ^ At{x)\XP-''+^ is in 
C''([0,T] X ^7^'-^S(X^'-'■+^X^'-'■)). 

Suppose there exists a linear subspace K of X, contained and dense in 
XP for all p > 0, and suppose that for all p > there exist G 
C^i[0,T],B{XP+\XP)) and operators {Ei''\x)}^t,x)€lo,T]xU^+^ in X such 
that 

(viii) 5'^^^'* : XP~^^ — > is an isomorphism for each p, and an isometry 
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at t = 0; 

(ix) ^Sj^^x' = [At{x), S^^]x' + E^\x)S^^x' for all x' e K and all {t,x) G 
[0,r] X [/P+i; 

(x) [0,r] X C/P+i 3 Ei^\x) is 5(XP)-norm bounded. 
Suppose also that there exists an integer q>2 such that 

(xi) the spectrum of {Aq)~^ Fo(uP) G B{X) contains no integer in the inter- 
val [l,q]] 

(xii) \\{{AI)-^Fq){vP)\\b(x^) < g - 1 for all p : < p < g. 

Lastly suppose that : [0,T] x ^7 ^ B{X) satisfies 
(xih) A^{x) has range X for all G [0,r] x U; 

(xiv) there exists a > such that al — A^^ix) is dissipative on ^(t^x) foi^ 
all (t,x) G [0,r] X U. (i.e IK^? + (A - a)/)x||x(,,,) > A||x||x(,,,) VA > 0); 

(xv) vlj(x) leaves invariant for all {t,x) G [0, T] x U and all p > 0; 

(xvi) UP 3 X ^ A^^{x)\XP is in Lip{UP, B{XP)) for aU t G [0,r] and 
all p>0; 

(xvii) [0,r] xUP 3 {t,x) ^ ^°(x)|XP is in CP{[0,T] x UP,B{XP)) for 
all p > 0; 

(xviii) [0,T] X C/P+i 9 (t,x) ^ [{Al{x))-^ , s[^\At{s[^'^)-^ is 5(XP)-norm 
bounded for all p : < p < — 1. 

Let 

C/O = {(^,0)/ e u , Ft{x){u^)' = V(t, x) G [0, T] x U} (14) 

Then there exists Tq G (0, T] and an open neighbourhood Jj'^i^^ c [7^'?+-'^ of n' 
in such that, for each (u^)' € f] U'^i+^ the evolution problem (13) 

has a unique solution u' G np=o ^^([0' ^o]) t^^~^)) and the mapping n 
^2q+i 3 (^oy ^ ^/ g (;70Qo, To], [/) is continuous with respect to the X^^+i 
topology on C/° n 

By taking XP = HP+%W\M.^) , where F''(M™, M^) is the L^-type Sobolev 
space (Adams 1975) of R-'^-valued functions on W^, and s > m/2, Newman- 
Claudel arc able to apply the above theorem to symmetric hyperbolic PDE, 
yielding the following result: 

Theorem 2 (Newman-Claudel 1998). Consider the system of quasi- 
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linear PDE 



^0 



{ut)dtut = A^{ut)djUt + f{ut) + -F{ut)uu 0<t<T 



n(0) = n' 



,0 



(15) 



on M™' for an unknown function u : [0, T] i-^ U, where U is an open neigh- 
bourhood of u*^ in X = H^{W^,R^), f is an X- valued function on U, 
and A^,... , A"^, F are S(X)-valued functions on U. 



(1) A^{x), . . . , A^{x),F{x) have the form oi N x N matrices, and f{x) has 
the form of an AT-component column vector with each matrix or column ele- 
ment a rational function (over M) of the column elements 
Xl, . . . ,xn of X] 

(2) A^{x),... , A^lx) are symmetric matrices. 

Assume also that the initial data satisfies 

(3) F{x)u° = for all xeU; 

(4) the range of n° in is uniformly bounded away from the zero sets of 
the denominators in (1); 

(5) A^{u^) is uniformly positive definite as an N x N real matrix- valued 
function on M"*; 



Finally assume there exists an integer q >2 such that 

(7) for all ^ G M"* the matrices (A°)~^F('U°(^)) have no integer eigen- 
values in the interval [l,q]; 



(8) \\{A^)-'F{u%0)\\b(r^) < 9 - 1 holds for ^ G 

Then for initial data u° G iJ29+i+s(M"^, R^), there exists a unique solu- 
tion u G nl^oCP{[0,To],Hi-P+'{W^,R^)) of (15), with an 
i7^'?+^+''(M™, IR^)-small change in the initial data giving rise to a 
C°([0,ro],if^(M"*,M^))-smaU change in u. 



Assume 




8 



4 The Cauchy problem for polytropic perfect fluid 
spacetimes 

In this section we consider the conformal Einstein equations for polytropic 
perfect fluid spacetimes which admit an isotropic singularity, extending the 
work of Tod (1990, 1991) and Newman (1993a,b). The hard work is to find 
suitable variables to put a reduced form of the equations into the form (15), 
and to see that a solution of the reduced equations does give a solution of the 
Einstein equations. The conclusion will be that for 1 < 7 < 2 the conformal 
Cauchy problem for the conformal Einstein equations is well posed near E, 
with the free data consisting of just the 3-metric there. The existence and 
uniqueness results obtained are used to prove the Vanishing Weyl Tensor 
conjecture described in the Introduction. 

4.1 Conformal gauge fixing 

Suppose then that Q is the conformal factor relating the physical metric gab 
to the unphysical metric gab, by means of 

gab = ^^gab (16) 

with ri = at S in M. Suppose also that the stress-energy-momentum 
tensor in M is that appropriate to a polytropic perfect fluid: 

Tab = {P+P) UaUb - Pgab ( 1 7) 

p={j-l)p 1 < 7 < 2 (18) 

To make analytical progress we assume that the rescaled fluid flow given 
by Ua = ^^^Ua remains regular and non-zero near E. 

From the work of Newman we know that if one demands that Q be smooth 
with non-zero gradient at E then the Einstein equations imply that the value 
of the polytropic index must be 7 = 4/3, the value for radiation. So if one 
wants to consider other equations of state then a different differentiability 
conditon on is required. To see what this should be, consider the FRW 
cosmologies, where one may take the scale factor R(t) for O. With a perfect 
fluid source, the field equations imply that the behaviour of dR at i? = is 
tied to the value of the polytropic index 7 according to 

dR ~ as r ^ (19) 
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in terms of conformal time r which is a smooth coordinate in M. In view 
of this, we will impose the following condition: 

At S , Z = n (37-2)/2 jg smooth, and VaZ ^ (20) 

It is a result of Goodc (1987) that the physical velocity field Ua must 
be hypersurface orthogonal, and this will also be true for the unphysical 
velocity field Ua = ^~^Ua ■ We know also that Ua meets S orthogonally. 
Thus near S in M 

Ua = ^Z,a = g"-''Z,J,b (21) 

for some smooth cosmic time function Z on M, with Z = onT,. A natural 
choice of conformal factor is then given by the following lemma: 

Lemma 3.1 If (20) holds then one may take the conformal factor Q to 
be a function of the cosmic time defined by the fluid velocity. 

Proof We know Ua = ip'VaZ for some smooth tjj, Z in A4 near S, 
with Z = at S. Now gab = Qah , where Qab = (U"^ /^'^)gab, and 
n = (Z)(i/5) = (3^ _ 2)/2) . So gab = {Z/ Z)(^/i) g^b , where Z is as 
in (20). Since Va-Z^, Va-^ are non-zero in M , one has Z = Z f , for some 
smooth function / in M, which does not vanish at S. Therefore fC^/i^ is 
smooth at S, and gab is smooth in M □ 



4.2 Conformal field equations in M 

We now make explicit the field equations we wish to solve for g^b in M. 
If we make the conformal gauge choice of Lemma 3.1, with Z as in (20) 

and Ua, V as in the unhatted version of (21), then for a physical stress 
tensor of the form (17) the conformally transformed Einstein-perfect fluid 
equations for g^b in M are (Tod 1991): 

4 1 127 1 

- §73^ + (37-2)2Z^^"^^^^ " 



Here Za = VaZ and G parametrizes the density according to 

12 V^G 

P~ (37 - 2)2 ^67/(37-2) 



(23) 
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This is the field equation in M for which we wish to formulate the Cauchy 
problem. 

Firstly there is information in the contracted Bianchi identity applied to 
(22). Projecting along and perpendicular to Ua gives, respectively 

□^^(2^£!>4^iz^^„ P4) 

7 V 7 G 

(2 - 7)^ - (7 - 1)§ + ((7 - 1)^ - (2 - 7)^) || = (25) 

where K = VaF, G„ = VaG. 

If 7 7^ 1 then (25) can be solved to give 

G = VPgiZ) P = (2-7)/(7-l) (26) 

for some function g. 

The case 7 = 1 will be treated shortly. Now we aim to show that g can 
be set to 1 by a change of cosmic time Z. Multiplying (22) by Z"^ and 
letting Z ^ we get G ^ 1. Also, we see that Z^^{1 — G) has a limit 
as Z — > 0. Now multiply (22) by ZZ^ and let Z — >^ to find that Vb is propor- 
tional to Zjj in the limit. This means that V must be constant on E and we 
can take this constant to be 1 wlog. From (26) this entails ^ 1 as Z — > 0. 

Suppose now that we do (16) two different ways: 

9ab = ^'^9ab = ^"^gab (27) 

then from (23) since these both correspond to the same p, 

^2^^^^67/(37-2) ^ y2^_^67/(37-2) ^28) 

Now use (26) to get 



^.(7-l)/7 Y /^N 6(7-l)/(37-2) 



V[-z) 



(29) 



Finally V/V = {n/n){dZ/dZ) so that 



dZ \g 



(30) 



Solving (30) for Z with g = 1 determines a new choice of unphysical 
spacetime in which g = 1 , and (26) becomes 

G = F(2-7)/(7-1) (31) 
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while (24) becomes 

UZ+ P-^ = (32) 

If we let hab = Qab — UaUh be the metric on the surfaces of constant Z , 
and Kab = \Cuhab the second fundamental form, then we have the identity 

where K = h^^Kab- 

Therefore (32) can be written as 

dV 
QZ= 

The field equations to be solved are now (22), (31) and (34) 
Returning to the case 7 = 1 we find from (25) that F is a function 
only of Z. Now we can use a change of cosmic time to set F = 1 . To 
see this note that if Z, Z are two choices of cosmic time with corresponding 
functions fl, V and tt, V then 



dZ _VQ _V f zV 

dz ~ vn~ V [z I 



(35) 



Prom (22) wc again deduce that V — > constant on S and we can take this 
constant to be 1. Now we may put V = 1 in (30) and solve for Z to get the 
required result. 

In this case (24) becomes 



and thus by (33) one gets 



UZ + — ^ = (36) 



and the field equations are (22) and (37) 

4.3 Solving 'reduced' field equations for 1 < 7 < 2 

If one wishes to give data for (22), (34) at S then the equations impose 
constraints. We've seen that V must be equal to one, and one also gets = 
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from (22) (or see (41) below). The equations impose no constraints on the 
initial 3-metric hab , and this is all the free data at S. 

The aim now is to cast the equations (22) , (34) in a symmetric hyperbolic 
form for which Theorem 2 applies. 

We will use coordinates Z, , with the being carried along the 
integral curves of V"Z, on a manifold of the form E x [0, T]. 

In these coordinates the metric g assumes the form 

ds'^ = ^dZ"^ + hijdx'dx^ (38) 

with hij negative definite. 

Making a 3 + 1 split of equations (22) gives 

= ^(4 - 37)/(37 - 2) + A(yP _ - 2) (39) 

i^. = |(^a.logy (40) 

P - 4 V ( 2V{2-^) K 6(2-7) p \ 

{3j-2)Z^'' + (37-2) Z + (37-2)2^2^^ ^ >) 

(41) 

The terms which appear here are undesirable with regard to ap- 
plication of the Newman-Claudel theorem. We therefore introduce a new 
variable C, according to 

6(2-7) {V^-l) 
^ " (37 - 2)2 Z ^ ' 



Prom (34) there follows 



^ = - 1 6(2-7)^ 

dZ Z^ Z (37 - 2)2 



(43) 



By definition one has 



Rab = + d{a^V) + ^ah (44) 

where 
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(Note that Hab is not the magnetic part of the Weyl tensor.) For a metric 
of the form (38) one therefore has 

Ti = di log V + /i™"(9^/lni - \dih„,n) (46) 

T, = -d,\ogV + ^hi,k'^ (47) 

where k^^ = (^ Kij = -^Vhi^hjnk'^'' , K = -\V{hijk'^) ). 
Also 



{di\ogV){dj\ogV) -r^.^fciogy + \v\d,\ogV)hkihjik 



kl 



Hiz = {hk\jdi\ \ogV)k 



jk 



1 



+ ^hijhmnihV^ - ^"^^F«)r^, (49) 



^(&iogy)/i™„,fc™" - -hi,hMk'''k^^ 



H,,= 2(0, log y)^ + -id,logV)hmnk"'^ - -hi.hkit 

+ ^h'^ {{di log V){dj log V) - V%dk log V) (50) 

Define ^ij^ = dj-hij . Then from the definitions (46), (47) one gets 

d.Ti = diV, - hikdj^^ + 2did, log V - hikh^^'^y^^mjn (51) 
Now (34), (47) imply 

= ^ (2 - l)hijk'i = Pd, log V (52) 

whereby (51) becomes 

dzTi = (^) ^^r. - hikdjk^^ - hikh^^y'^^mjn (53) 

We now choose new independent variables in such a way that the field equa- 
tions for these quantities become symmetric hyperbolic. 
Define the positive constants E, F, G by 

2 _ 247(2 -7) p2_ 127 ^2 _ ^7 
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and define the following new variables 

Xi= E(di\ogV + ^T^ , v= FT, , ii = GVi (55) 



In terms of the new variables the field equations can be written as follows: 

d,W^ = (56) 
dzhij = h'lmhjjik (5'^) 

^ 2(1 - 57) V ^ ' (2 - 7) i^y 

z (37 - 2)2(1 - 57) V (2 - 7) i^y y ^ ^ 

+ G /i"'"/i;'^7mZn + 2G/l*^iIi, 

+ 2V'^k''"^k^h(^i^)hQp)h(^kr) — "^Hij 

/ 2 / /''(mn)^™" 2(7-3) 



Z V (87-2) ^""^ ^^"^ -w; (37-2) ' (37-2)(7-2)F 

(62) 

(2^)5..= ^Ig^'^'^^^^- 

+ 2F/l*^' -M-] + _ ^2^^ 



£; 2PG ME' 2PG / FP2 
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where 



/ /37_2\ h,,„„J."- 2(3,-5) .\ 
+ Z VV2-7r+ (37-2) +(2-7)(37-2)W ^''^ 

^ (37-2)Pz'' \2PG e) ' 



+ hijhmnihP'^k"'' -h"'^kPi)t;g (66) 

Xi _ (Xl_ \ _ f fe /^Xfe _ 

E 2PG J \ E 2PG J '^\E 2PG 

+ lv^hk^ihj)ik''"'k'''hmn+ \v\j-2)hmnk"'^hk^Jl,)ik''' (67) 

rf, = f f,^.) ; f ^. = h>'\ju, + 7,7. - 7^,7) (68) 

In (58), (59), (62), (63) we have used the first of (52) to partially sub- 
stitute for hijk^^ , which will be useful later on. Equation (62) comes from 
(44) and (41). Equation (63) comes from (44), (39), and the second of (52). 
Equation (61) comes from (40), (44), and (53) was used to partially sub- 
stitute for diTz- Equation (64) comes from (44), (40), and the second of (52) 
We now seek a solution of the 'reduced' equations (56)-(68), 
where 'jijk, Xi, i^, ^^e to be treated as independent variables, with 
no reference to their original definitions. We will suppose that the coor- 
dinates are chosen to satisfy the harmonic gauge condition 

h'^r'^j = (69) 



and 
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on S. 

The initial data for the reduced equations are as fohows: 
hij is free data at S , and is chosen negative definite and symmetric, h^^ is 
chosen so that h^^hjk = and choose also 7jjfc = d^hij. Since the second 
fundamental form vanishes at S one must have k^^ = there, and one 
also has y = 1. Prom the first of (52) and (53), (56), (55), (69) there 

follows V = ( = Xi = ?i = 

It is possible, using the formalism of Newman (1993b), to write (56)-(68) 
explicitly in the following form 

A°{u)d^u= A\u)di{u)+ B{u)u+ ^C(tx)u (70) 

where u stands for the field variables written as a vector, is positive 
definite and symmetric, the A!" are symmetric, and A", B, C are polynomial 
in their arguments (see Appendix for details). One has C(m)(u(0)) = for 
all u, and the eigenvalues of ((^'^)~^C)(«(0)) satisfy one of the following 

A = 0, A- 



P(37 - 2)2 

_^3^ 3(2-7) ^^2 , 12(2-7)(57-3) ^ 32(137-8) 



(3^ _ 2) (37 - 2)2(57 - 1)) (57 - 1)(37 - 2)3 

Hence for 1 < 7 < 2, no A is positive. It follows that (56)-(68) satisfy all 
the requirements of Theorem 2, and thus there exists a solution u unique in 
a suitable differentiability class. Differentiability is discussed at the end of 
section 4.4. 

One now wishes to know whether the metric g obtained from u via (38) is 
actually a solution of the conformal Einstein equations (22). i.e. whether 
the definitions of 7yA:, Xii ^) can be recovered from (56)-(58) . The 
answer is yes, as is shown in the next section. 



4.4 Recovering the conformal Einstein equations 

We now show that if u is the solution of (56)-(58) as above, then as long 
as remains negative definite, is in fact symmetric, and g given by 
(38) is a solution of the conformal Einstein-perfect fluid equations (22). The 
key step will be, following Newman (1993b), to use the contracted Bianchi 
identities to show that vjF and ^ijG agree with the definitions (46)-(47) 
of r^; and Fj. 
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Suppose then that hij is negative definite. Equation (62) then gives 

a,(Z4/(3T-2)fc[ii]) = 

Hence fel*^] = for all Z. Then (56)-(57) give hiij] = h^'^^ = for all Z. One 
also gets 

d.ih^nij^ - 5f) = -hUh^'^hjn - St)k^^ 

A Gronwall estimate (Racke 1992) therefore gives h'^^hji = 5^ for all Z. Now 
(56), (60) imply 

dzilijk ~ dk^ij) — hjmk™'"' {''fink ~ dkhin) + hirak^^i^njk ~ d^hnj) 

So a Gronwall estimate implies ^yijk = dkhij for all Z. 
Define ^ by 

6(2-7) jV-l) 
^ (37 - 2)2 Z 

then (58), (59) give 

d,{z{c-0) = o 

and hence the definition (52) is recovered. 

Now write t^, = ^ , fi = §, V=§- Then (61), (64) imply 



2-1 J ' \ 2t ) • (7-1) 



and now (58) gives 



^ + 2PVi + f ^ir^ ) - 2Pdi log y = (71) 



27 y " V 27 

Note that once we have established Fj = Fj then (71) will imply 
V'i = 5i log y + 2^rj as desired. 

Prom the definition of and (58) there follows 
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while (71) implies 



- 2^r, = (r, -v,) + di log V (73) 

and (64) implies 

_a,,. = -A^8,f, + i4. + M5^(2^f.-*) ,74) 
Equations (72)- (74) and (51) now combine to yield 

Define a symmetric tensor Sab by 

~ Z \ (37-2) ^'"^'^ 



4 1/ 1 

V'^-TTT^rJ (77) 



Z (37 - 2) V 2P 

z\v2-7y^ 37-2 (2- 



2(37 - 5) 



7)(37-2) 



(78) 



Also define Tab = Sab - \Sgab , where S = g^'^Sab- 
Then 

_ 1 r 2 4 .1 

^^^-zl(2^+(2-7)(37-2)^^/ ^^^^ 
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Equations (62), (64), (63) can now be written, respectively, as 

Sij = Rij + %(rj) - r^)) + {Hij - Hij) (82) 

Siz = Riz "I" {Hiz H-iz^ 

+ l (l^) - r.) + (^^) dziti - Ti) (83) 



Szz = Rzz + {Hzz — Hzz) + 2^3^ _ ^z(yz — r^) 

+ (\^) /^-(r.-f.)a,iogy+i (3^)^r.-r.)^+ ^^^(^^-/2;^^ f.(r.-f.) 

(84) 

If (56)-(68) hold, then one calculates the following 

+ (^^) a >"8 ^'"(f'- - Tz) + (^) (r, - fi)hirak'"' 

, 3 (2 - - i)(r.-f.)+(2fi) f f i±l2) r, + i'-^'^jl f,) (r.-f.) 



Z 2(37-2) ' ' V 47 7 V\l-37y (2-7)(l-37) 

(85) 



V T- = Yl\ \ / 2(137-8)-6(y^-l)(57^-47 + 3) \ 

^ (l-37)(37-2)2 ; 

+2f (f D (^-^)(^^-^) ] 

+2i,(i, i^)(3^_2)(2-^)(i-37)/ 

4 { (i^^'^(^^- - + (3^ icw)^^'^^' - ^'^^^^ - 

(86) 
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Now write Aab = Sab - Rab = {Tab - \Tgab) - Rab- Then by the Bianchi 
identities one gets 



/*9,A„5 = V,Tl - \da{T + R) + 2/''rJ,A,), (87) 



2 

and note that T + R= -g^'^Aab- 
Prom (82)- (84) A„b is given by 

Aij = d^iitj) - r^)) + terms in (f „ - r„) (88) 

= I iw^) ^^(^^ - + (^) ^^(^^ - 

+ terms in (f„ - r„) (89) 

+ terms in (f„ - r„) (90) 
Now substituting (85), (86), (88)- (90) into (87) gives 

= ¥ (2^) ' - n (^) ^-^^ff » - ^») 

- ^'"'I'Z w~ ^' ''^'^Z'-' + teriM in(f„-r„). z-i(r.-rj, aj(f„-rj 



/^5cA,, = ^ a.2(f . - r.) + (i±^ ) h--dUrn - r.) 

^^\r _p ^ ^2(137-8) -6(F^-1)(572- 47 + 3) ^^ 



^ (i_3^)(3^_2)2 ; 

+ terms in (f„ - r„), Z~\ta - TJ, ^^(f „ - (92) 
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The unspecified terms in these equations are polynomial in the stated quan- 
tities, with coefficients which are polynomial in the components of u. 

One can also calculate g'^'^dcAab directly from (88)-(90), giving 



+i/i-9L(f ^ - r.) + \h^-dUtr. - r„) 

+ terms as in (91) (93) 



/^5eA,, = ^ (^^) dl{v, - r,) + j ^"^"5L(f n - r„) 

+ terms as in (92) (94) 
Now we can compare (91) with (93) and (92) with (94) to obtain 



2 ' ' Z 1 (37-2)2 I Z 

^^((|^)«.tf.-r...(^)a.ft-r. _ 

+ terms as in (91) (95) 



dUt, - r,) = %^/.--5L(fn - r„) 



4 V37-iy 87 
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1/1 + 27 



2 V37- 1 



^"""^mnlr^ - r^) + terms as in (93) 



y2 r^-VA / 2(137-8)-6(y^-l)(57^-47 + 3) \ 
+ -Z ) [ (l-37)(37-2)2 ) (^6) 

Next we introduce new variables as follows 

■Ki = ti-ri, T = tz-Tz, aij = di{tj -Tj), Pi = dz{ti-Ti) 

where 



27 + 1 / 27 



2(37 - 1) ' V 37 - 1 

In terms of these new variables (95) becomes 



47 y ^"^^ 4(37 -i)C2 V 87 

+^h^-' h!'^ dkami + order zero terms 

V\,^\ 3(2 -7) (7-1) 2(3-47) (3-47) ^ 

Z'' \ (37-2)2 ''^^ (37-2)(l-37)Ci *^7(37-2)^' 

(97) 

and (96) becomes 

T (^) - ^'-^ (^) ~ - ^ (^) — 

^l^^ 2(137 - 8) - 6(y^ - 1)(572 -47 + 3)^ 



Z V (l-37)(37-2)2 

+ order zero terms (98) 

By the above definitions one also obtains the following 

d,-Ki= A (99) 
d,T= [C^r^b (100) 
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(101) 

-V^h'^d,Uj= -—^h'^djS (102) 

-V\'^d,rjj= -Ylh'^{-r,j + Pj) (103) 

V%K= }^{{C2)-'S - k) (104) 
The system (97)- (104) is of the form (see Appendix): 

a'^{v)d,v= a\v)div+ b{v)v + ^c{v)v (105) 

where v stands for the variables written as a vector, is symmetric and 
positive definite, the a* are symmetric, and aU coefficients are polynomial in 
the components of v and those of u ( u as in (70)) 
The eigenvalues of {aP)~^c{0) satisfy one of 

(l + AjfA + fl^V 12,(2 -,)(,-!)_„ 



, (37-2)V (37-2)3 
Hence for 1 < 7 < 2, no A is a positive integer. 

The initial data for (97)-(104) is as follows: One has Fj = at Z = 0. 
Since cp is constant at Z = one also has Tj = there. Hence tTj = aij = 
initially. Now = z/ = at Z = 0. Since % = at Z = 0, (58) 
implies dzV = there. Hence = at Z = and so r = Ui = there. 
Taking a Taylor expansion of (61), (64) gives dzS,i = dzTi = at Z = 0. 
The identity (51) then gives S^Fj = at Z = 0, and so A = there. 
A Taylor expansion of (59), (63), (58) gives = ^d^tz at Z = . 

Comparison with (72) then gives dz{Tz — ^z) = at Z = 0. Thus S = k = 
at Z = 0. 

All the initial data for (97)-(104) are therefore zero, and one solution 
of this system is w = 0. For sufficiently smooth v this solution is unique 
by Theorem 2, and we have recovered F^ = (An analogue of Theorem 2 
can be obtained from Theorem 1 in the case when certain of the coefficients 
are polynomial in suitably regular known functions) . The metric g obtained 
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from u in (70) via (38) is thus a solution of the conformal Einstein-perfect 
fluid equations (22). 

We now discuss the differentiability of solutions to (56)- (68) and (97)- 
(104). Suppose then that S is a smooth paracompact manifold. In or- 
der to apply Theorem 1 one must first consider the localised problem, and 
then use the finite speed of propagation inherent in (70), (105) to glob- 
alise (see Claudel and Newman 1998). Specifically, choose a locally finite 
open cover {Oa} of harmonic coordinate charts for the metric hij on S, 
with each Oq. having compact closure therein. And let {O'j} be another 
cover for S satisfying O'^ C Oq. Suppose that the initial metric lies in the 
following Sobolev class 

{hij)ae , s>3/2 (106) 

where 

ga >max{3+||(aO)-^c||s(M24), l + \\{AY'F{0,x)\\B(Re3)} , x e Oa (107) 

Then Theorem 2 gives that there is a solution (u)'^ of (56)- (58) contained 
and unique in the following class 

«e r]C\[0,T,J, m'^+^-\0'^)) (108) 

1=0 

Now by the definition of v in (105) there follows 

«G n C\%T,^],m'^+^-^-\0'^)) (109) 
1=0 

and (an analogue of) Theorem 2 gives that v = . Hence the met- 
ric {gabYa obtained from (u'^) is a solution of the conformal Einstein-perfect 
fluid equations (22), and one has 

{gabYa & nC^([0,r,J,if««+-^(O;)) (110) 
1=0 

with {gabYa unique in this class. 

Note that if we give C°° data {hij)a at S then the evolution given by (108) is 
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C°° on some time interval [0,Ta]. For the data appropriate to the locahsed 

problem is Cg° and thus one gets (u)^ as in (108) for arbitrarily large q^. 
For fixed a the Tq^ are positively bounded below by the following standard 
extension argument: first from (108) one gets by the Sobolev imbedding 
theorem (Adams 1975) that {u)'^ is on [0, T^,] for some > 0. Suppose 
now that (n)^ € C{[0,Ti), H'') for some r > |, Ti < r„. Then, by propo- 
sition 5.1.E of (Taylor 1991), 11(^^)^11//'' is bounded on [0,Ti) and (u)'^ can 
be uniquely extended in onto [0, Ti]. Then standard theory for regular 
symmetric hyperbolic PDE gives an extension of (u)'^ onto a larger inter- 
val [0,r2) with 

r 

{ny,ef]C\[T,,T2),H^-\0'J) 

1=0 

Thus (u)'^ can be extended onto [0, T^] in the above way. It follows by the 
Sobolev imbedding theorem that (u)'^ is C°° on [0, Tq]- 



5 Solving the Cauchy problem for 7 = 2 

The method of section (4) is clearly inadequate for dealing with the cases 
7 = 1 and 7 = 2 as certain of the coefficients in the field equations diverge 
at these values of the polytropic index. However, in the case 7 = 2 one 
must have HZ = by (26), (32), and it turns out that the problem may 
be successfully attacked using harmonic coordinates in M. Unfortunately 
in the case 7 = 1 it has not been possible to write the field equations in 
symmetric hyperbolic form, and the problem has not been solved in such 
generality. However in section 6 we are able to solve the problem in the case 
of Bianchi type spatial homogeneity. 

For the case 7 = 2 first note that the conformal field equations are just 

VaVbZ = ZRab (111) 

Again we wish to solve the Cauchy problem on M = E x [0, T] . We will use 
coordinates = Z, a;*, i=l,2,3 on M, and in the sequel Greek indices will 
take the values 0, 1, 2, 3. 

Now let iJ" = nx°', and consider the following reduced equations for g^i, 

^<. = V^V.Z--^A^, (112) 
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where 

R^, = R^.u+ ga(^d,)H^ (113) 

Aij = Aio = Aoi = 0, Aoo = 1 (114) 

Note that (112) impUes = so that Gjj^, = R^^ 

Now write hap^ = d^gap- Then (112) can be expressed as follows 

- g^''^d,hapj = -g^'^^djh^af3)0 (115) 
dzga/i = ha/30 (116) 

5.5"^ = -^"'^/'^Vo (117) 

g'^d.h^po = -g''''^dih^ap)j - 2g''d,h^a/3)o 

+ 2-ff(a/3) + -^g^'^i-ha/Sa + hpaa + haa/s) (US) 

(for a, P not both zero) 



1 . 1 ai.,. . 2 

00 + 

where 



+2Hoo+^{^g''g''\hijk-2hikj)+^g''''g''\hoki-hoik-hkio)+Q} (119) 



Q = g'Hhijo ' 2h,oj) - 25°'/ioo» (120) 

We can write the singular part of (118) as ^Ci{u)u in such a way that the 
components of Ci{u)u are 

g^^{h/3Qa + haQ/3 — ^a/3o) + (^/3iQ + " ^a/3i)5°* (121) 

where the tilded terms are in Ci(n) and the untildcd terms are in u. And 
the singular part of (119) can be written as -2C2{u)u with components 

1 ~ ~ . . 2 ~ 

-00 (^y^ ~ ^hikj)g^^g + -^{hoki — hoik — hkio)g ^g 

+ g'^ ih,,o - 2hio,) - 2/'/?,oo« (122) 

where the tilded terms arc in C2(m) and the untilded terms are in u. 
Clearly equations (115)- (119) are of the form: 

A^{u)dzu= A\u)di{u) + B{u)u + ^C{u)u (123) 

where u stands for the field variables written as a vector, the are symmet- 
ric with positive definite, and the components of C are those of Ci, C2 writ- 
ten in an appropriate order. 
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5.1 Initial data for the reduced equations 

For the initial data we first choose a negative definite 3-metric hab on S. 
Next choose coordinates on S satisfying the harmonic gauge condi- 
tion hP'^r^jf, = 0, and also take gQi = 0. We know from the field equations 
that g^^ = 1, hijQ = initially, and hence 500 = 1 initially. We make the 
choices /laoo(O) = 0, and this ensures Dx" = at S. These choices also 
ensure that {C{u))u{0) = Vn. 

5.2 The eigenvalues of (A°)-^(7(0) 

The eigenvalue equation for (^°)~^C(0) amounts to the following 

g^^(hi3ia + haijd — ha/3i) 

(with hooo replaced by g^-'{hijo — 2/ijoj) in rhs) 

^hapk = 

Xg"' = 

(Unbarred quantities are elements of an eigenvector, barred quantities are 
evaluated at S ) 

So suppose A 7^ 0. Then 1 implies 

^hijo = hijQ + g {hjmi + himj hijm\ 

Since the are harmonic in S this implies 

^'St^hijQ = -g'^hijo 

Prom 1 one also gets 

-^^000 = 9^'' hijQ , A/jjoo = 

It follows that A = or —1. The system (115)-(119) therefore satisfies all 
the requirements of Theorem 2, and has a solution u for each choice of 
hah on S, unique in a suitable differentiability class. Our next aim is to 
show, as in section (19), that the 4-metric gab obtained from u is actually a 
solution of the full field equations (111). 
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5.3 Recovering the Einstein equations 

Firstly we show that if (115)-(119) are satisfied, then g^is is symmetric and 
the definition ha/3j = djga/3 can be recovered. 
Observe that (118)-(119) imply 

and (115) implies 

Thus h[ai3]j = ) aiid there follows 

dz{Zhio,i3]Q) = 

Hence /i[„a]o = , and now (2.101) implies gai3 = g/Sa- 
Prom (115) one now gets 

dz{ha/3j - djga/s) = djha/30 - dljgafi = 

by (116), so that h^pj = djg^p. 

We will now use the Bianchi identities to show that if equations (115)- 
(119) (or equivalently (112)) are satisfied then we have Da;" = , and thus 
(111) is solved. 

So suppose (115)-(119) hold. Then one calculates 

+ g^^ {^,,d,{H'/g'') - {H^/g''){Tl^Ko + r°.Ao^)} ] (124) 
The Bianchi identities now give 

= g^^dlH"^ + 2g^'dliH'^ + g'^d^H'^ 

+^[g'^'d,H- - g'^-d.H' + 2g'-g'^{{g'^)-^d,H^ - {g'')-'H%oop) 

+ 2(ff0"/5°°)(i?°)2 - 2<7'^-ffO^(i?0/<7°°)r°,} (125) 

If we now put ha = dali^ , then (125) can be written in first order form 
as follows 

- g'^d.h^ = -g'^djK (126) 
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d,H^ = (127) 

- 2(/«/500)(ij0)2 + 25'^"/''(iJ°/5°'^)r°,} (128) 
This system is symmetric hyperbolic: 

a^{v)dzv = a\v)diV + b{v)v + ^c{v)v (129) 

where v stands for i?"^, h^. We know that iJ" = at S , and now equation 
(128) impHes that /iq = at E also. Hence all the initial data for (126)-(128) 
are zero, and one solution is trivial. 

The eigenvalues of c(0) satisfy either A = or 

A« = -h% + - 2f"fPhl 

Putting a = 0, j gives, respectively 

Xhl = -2/i[] , Xhi = -hi 

Hence A G {0, —1, —2}. So (126)-(128) is in Newman-Claudel form, and the 
zero solution is unique in a suitable differentiability class, so that (111) is 
satisfied. 

5.4 Existence, uniqueness, and differentiability 

The situation here is very similar to that in section (4). So take cov- 
ers {Oa}, {Oq} for E as before, and suppose that the initial 3-metric lies in 
the following class 

{hij)a e H^i'^+^+%Oa) , s > 3/2 (130) 

where 

qa> max{3+||(a°)-ic(0,a;)||s(R96) , l + \\{A°)-^C{0,y)\\B(R2o^} x,yeO^ 

(131) 

Then by Theorem 2 there exists a solution u of (112) contained and unique 
in the following class 

«G nC^([0,r,J,i7««+-^(O;)) (132) 

1=0 
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Prom the definition of H°', ha it follows that 

{?a-2) 

«G n C^X[0,r,J, i?««+-2-(0;)) (133) 

i=0 

and by (an analogue of) Theorem 2 one must have v = 0. Hence the metric 
gab obtained from u is the unique solution of the Einstein equations (111) 
with starting metric hij on S. 

We note that, as in section 4, C°° data gives rise to a C°° solution u. 
5.5 The Vanishing Weyl Tensor Conjecture 

The results of sections 4-5 can now be used to prove the Vanishing Weyl 
Tensor conjecture. From (Goode and Wainwright 1985), and as one can see 
from equation (9) in this case, one knows that the Weyl curvature vanishes 
at S iff the initial 3-metric is of constant curvature, and that any such 
metric is initial data for some FRW cosmology. The 4-metric of any FRW 
cosmology is a smooth evolution of smooth data, and hence by the above is 
the only smooth evolution of the same data. Thus in this case the conjec- 
ture is true: for 1 < 7 < 2 at least the only cosmologies with an isotropic 
singularity and vanishing initial Weyl tensor are the FRW models where the 
Weyl tensor vanishes everywhere. 



6 Spatially homogeneous 7 = 1 spacetimes 

We noted above that for 7=1 the conformal Cauchy problem has not been 
solved in full generality. However if one restricts attention to certain (Bianchi 
type) spatially homogeneous cosmologies, the field equations become ODE's 
and a theorem of Rendall and Schmidt (1991) can be used to obtain the 
required existence and uniqueness result, as we shall see in this section. 

6.1 Bianchi type spacetimes 

The Bianchi type spacetimes are a class of spatially homogeneous cosmo- 
logical models having a 3-parameter Lie group G of isometries. They are 
classified into 9 types according to the nature of the Lie algebra associated 
with G (see Wald 1984 for details). In such spacetimes there exists a cosmic 
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time function t, and covector fields (e*)a such that the metric can be written 



as 



9ai. = ^atVbt + hij{t){e%{e^)b (134) 



The (e')(j are preserved under the action of G, and are called left-invariant one- 
form fields. 

6.2 Conformal field equations in M 

We now consider equations (22), (37) in the case where the physical metric 
gab has Bianchi type symmetry. We may assume wlog that the singularity 
surface is at Bianchi time t = 0. One must also have that the fluid velocity 
is orthogonal to the homogeneous hypersurfaces, by the following argument: 
There exist vector fields (^j)'* {i = 1,2,3) on the Lie group G which are 
Killing vector fields for gab{t), t > 0. Since the fluid velocity is geodesic 
for 7 = 1 there follows Uai^i)"" = const, in M. But Ua = , with Ua reg- 
ular. Letting ^ we see that Uai^i)"' = in M as required. 

The conformal factor is 0, = and the unphysical metric can be written 
gab = VaZVbZ + hij{Z){e')a{e% (135) 
since V = 1. It follows that $7 = (3t)^/^ where t is the Bianchi time. 

Suppose now that we write equations (22), in an obvious way, as Gab = 
Tab- In order that these equations be satisfied at Z = one must have G(0) = 
1 and KijiQ) = where Kij = ^dzhij. The initial 3-metric hij{0) is free 
data, as for the other polytropes. 

The equations Gij = Tij take the form 

dzKij = '^Rij - KKij + 2KimK;p 

-±K,-H,{^,%^) (136) 

d,h,j = 2Kij (137) 
= h^^Kij and ^Rij is tl 

given by 



where K = h^^Kij and ^Rij is the Ricci tensor of the surfaces Z = constant. 
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- \c2CLh,rnh,rh'''h'' (138) 

where Cj^ are the structure constants of the Lie group G. 
Now write 

C = {Gab - Tab)Z''Z' (139) 

and 

Cd = {Gab-TabW^Z' (140) 

A rather lengthy calculation then shows that if one has a smooth solution 
hij{Z) of (136)-(137), then the following hold: 

4' 



c{k+-) (141) 



d,Ci = -CAK + ^] (142) 



z, 

where Di is the derivative operator associated with hij. 
Now C and Q vanish at Z = by the conditions already imposed. If we 
regard K{Z) as a known smooth function then Theorem 1 of (Rendall and 
Schmidt 1991) implies that (142) has a unique solution d. Hence Ci is 
identically zero, and then in the same way C is identically zero. Our task is 
therefore to solve equations (37) and (136)-(137). 

It will be convenient to work with hij,h^^ ,G,K--' as independent vari- 
ables, and we also introduce the quantity (" via 

C = 6(1 - G)Z-^ (143) 

The field equations to be solved are then 

d.hij = 2h^,\k\Kif (144) 

d^h'^ = -2/i(*l'=li^/) (145) 

d,K/ = h^\^IUk) - K^K,^ - - \5,\2K^ + (146) 

d,G = -GK^ (147) 

a,c = -ci^r + ^(6i^r-c) (148) 

Equations (144)- (148) can be written in matrix form as 

^ + }-Nu = G{u) (149) 
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where u = {hij,h'^^,G, K^■' ,Q, N is a constant matrix, and G is polynomial 
in the components of u. 
The matrix N takes the form 

where the O's stand for blocks of zeros of various sizes and M is a square 
matrix. 

The eigenvalues A of M satisfy the following 

- AK/- = - S/{2K^ + (151) 

-XC = 6K^-C (152) 
Taking the trace of (151) gives 

(-A + 10)K„- = -3C (153) 
and then from (152) there follows 

(-A + 10)(-A + 1) + 18 = 0, or C = (154) 

Thus all the eigenvalues of M have strictly positive real part. Now if all 
the eigenvalues of N had strictly positive real parts, then Theorem 1 of 
(Rendall and Schmidt 1991) would imply the existence of a unique smooth 
solution u to equations (149) with the given data. However, in the case 
where N takes the form (150), with M having eigenvalues with strictly +ve 
real part, a very elementary modification of their proof gives the same result. 
It follows that equations (149) have a unique smooth solution u. 
It remains to show that if (149) is satisfied then hij = hji, h'^-'hjk = 
S\, K,, = hjkK..^^ = \d,hij, and C = 6^-^(1 - G). 

First note that (144), (145) imply that d^h^ij^ = d-^h^'^"^ = 0, so that hij, K-^ 
stay symmetric if they start that way. Also (144)- (145) imply 

d,{h^ihju - S^k) = ih^'hjm - - {y^hrni - 5^^)K^ (155) 

Hence by Gronwall's inequality h^^hj^ = if hij{0), h^^{0) are chosen as 

matrix inverses. 

By (144), (146) one has 

d.K^ij] = -{K^ + 4Z-i)i^[,,.] (156) 
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so that by Theorem 1 of (Rendah and Schmidt 1991) = and now 

(144) implies that Kij = ^dzhij. 
Finally, by (147)-(148) there follows 

d,{C-0 = -{Z~' + K){C-0 (157) 

so that C = C- 

We conclude that for each Bianchi type 3-metric hij on G , there exists 
a uniqTic Bianchi type 7 = 1 cosmology having an isotropic singularity, 
with hij as initial metric in the chosen gauge. In particular, the Vanishing 
Weyl Tensor conjecture goes through as for 1 < 7 < 2. 

7 Concluding remarks 

The idea of constructing cosmological singularities by means of an infinite 
conformal transformation provides a mathematical framework for studying 
the Weyl Tensor Hypothesis of Penrose. In particular one may seek to 
generate cosmologies with isotropic singularities via a Cauchy problem with 
data given at the singularity surface. By allowing the conformal factor to be 
non-smooth at the singularity we have been able to solve this problem for 
perfect fluid models with equations of state other than that for radiation. 
For these other equations of state there are as many cosmologies as for 
7 = 4/3 and they arc determined by the intrinsic geometry of the singularity 
surface, with no free data for the matter. 

If one wishes to impose the condition Cj,^^ = initially then the unique- 
ness results obtained show that the spacetime geometry must be exactly 
FRW in a neighbourhood of the big-bang singularity. 

One may now ask whether the picture painted above remains the same 
if one considers other matter models. What, for example, happens to the 
conjecture regarding vanishing of the Weyl tensor if we are able to give some 
matter data at the singularity? Indeed, are there matter models for which 
it is possible to give data at the singularity? This question will be addressed 
in a second paper in which we shall consider the Cauchy problem for the 
conformal Einstein- Vlasov equations. 

A Explicit matrix representations for the 1 < 7 < 2 
field equations 
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A.l Newman's formalism 



Here we present the formalism of Newman (1993b) which allows tensor ial 
evolution systems to be written in matrix form. 



Let V be an n-dimensional vector space with dual V*. Let V be equipped 
with a basis, and V* with a dual basis. Let h and h* be covariant and con- 
travariant tensors over V of rank 2. For any tensor T of total rank r over V, 
with components T'^ " ■ define the component row and column vectors 



.11) ■ 



.12) • • • ) ■ 



.In; ■ 



.21) • • • ) 



rpn. 



.12 

.In 
.21 



th 

A general matrix will be represented in the form M\, with M*. the i row 
and M'j the j^^ column.. For any x square matrix M'., and each i = 

1, 2, . . . , n let (*)M'. denote the matrices formed by the 1^^, 2^^*^, . . . , groups 
of n^^^ rows of M', respectively, and let (i)M ; be the matrices formed from 

the 1^*, 2^*^, . . . , n^^ groups of n^~^ columns of M',. Recall that the tensor 
product of a p X (/ matrix M\ and an r x s matrix N\ may be represented 
as a pr X qs matrix: 



( M\N\ M\N\ 
M\N\ M\N\ 



M\N-. \ 
M\N\ 



\ M\N\ M\N\ ■ ■ ■ M\N\ J 



Clearly this representation respects the associativity of tensor products. For 
each A; = 1, 2, 3, . . . one may therefore define the n'' x n'' matrix 



k 

ih)-. 



{hy.®...^{hy. 



k factors 
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where 



{hy. 



( hu hi2 

\ hnl hn2 



hln \ 
h2n 

hnn ) 



If is a symmetric tensor then [ji)' . is a symmetric matrix. Inductively one 
can show that if h and h* are tensor inverses in the sense of 

h*'^h^k = 5\ 

k k 

then (ji)' . and {h*)' . are matrix inverses: 



Y.^h*yM'. 



<T=1 



{h*y . is defined in terms of h* by analogy with the preceding definition 
k 

of {hy . in terms of h. Moreover, letting hg and /i* denote the symmetric 

k k 

parts of h and h* one may, for each A; = 1, 2, 3, . . . define {hg)' . and {h^y .. 
Suppose h and h* are tensor inverses and let T be any contravariant tensor 
of rank r over V . Then, defining a covariant tensor of rank r by 



hi 



hi i T-'i-^ 



^l\...lr 'Hljl 

one has a matrix representation of the raising and lowering operations: 



a=l a=l 

There exists a unique symmetric x matrix J', such that, for every 
second rank covariant tensor Tab i the tensor T^b = Tba satisfies 



(T)- = J■,(T)^ (f ) . = (r),j- 
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In particular, for n = 3 one has 



J\ = 



( 1 


































1 



































1 











1 



































1 



































1 











1 



































1 











V 























1 / 



J', also plays an analogous role for second rank contravariant tensors. 

r 

For each r = 1, 2, . . . there exists a unique rf x rf matrix E' . such that, for 
any contravariant or covariant tensor field T of rank r, one has 

(T,)- = e\ {Tf 

1 

where Ts is the symmetric part of T. In particular, for n = 3 one has E' = 
P and 



E = 



1 




























1 

2 





1 

2 























1 

2 











1 

2 











1 

2 





1 

2 





























1 





























1 

2 





1 

2 











1 

2 











1 

2 























1 

2 





1 

2 





























1 / 
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A. 2 The reduced field equations for 1 < 7 < 2 

If we define the 63-component column vector 



u 



( m- 




{hy 




V 




c 




w 




{ky 




(7)- 




(6- 




V ^ 


) 



then the 63 x 63 matrices A'^{u), C{u) appearing in the reduced equations 
(70) can be written 
A^{u) = 

/I' 



\ 



79 



o 



v^hsy 



o 



V 

A'{u) = 
/ 





(2-7)^ / 













-my- 








2 

E- 



.1 (Or- 
G ^ - 





{Ky.® E\ 







2 

G (i) 



2 



-msY 
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C{u) 

( 







where 



-/I 









)■ 




3(2-^)2 














D^[h%)- 




D^{\f:)- 




(37-2)2(1-57) 

_4P-1 



(37-2) 



Da 



(37-2)2(1-57) 



L>7 



(37-2) 

8GE-^ 



Da 



(2-37) 



p2(3^-2) 

(37-2)(7-2) 
4P-1 



(37-2) 



D.= 



y2F(37 - 2) „ F^F 
Dw = 



(2-7) 



(37-2) 



11 



2y2(3^ _ 5) 
(2-7)(37-2) 



A. 3 The recovery equations for 1 < 7 < 2 

Define the 24-component column vector 



/ (^)' \ 
(«)■ 

(^)- 
(5 



u = 



\ K J 



D2l^ 



DeiKY 



Di,I^ J 
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then the matrices a°'{v), c{v) appearing in the recovery equations (105) are 
as follows 



/I' 



2|'^,-l^ 



o 



o 



-v\h-^y. 



where 



while a*(i') 



7-1 ^^_ F^(2-37) 
87 47 

C5 = ^57 -r 



4 V37-1 



4 rl 



: 



QV(« ^-')" 












Ce{h-^y, 

C7(/i-i)-, 



: 
... 0/ 



where 



Q = -/9 + C3 J 

(3-7)^^ ^ l^V27 + AC^2 



4(87 -i)C2 2 Vi-37y Ci 
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and c{v) 
( 



where 





Cii/^ 

-v^{h-^y. v^{h-^y. o 

-VH^ J 



<^8 = -77; ^t; — = 



(37-2)7 Ci(37-2)(l-37) 



3y^(2-7)(7-l) 12(137 -8) + 6(1-^0(57 -47 + 3) 

^10 = 77, 7^ ^11 — ^ 



(37-2)2 - [ (i_3^)(3^_2)2 

and Ci, C2 were defined in section 4.4. 
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